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Abstract 

We present a one-parameter family of quantum maps whose spectral 
statistics are of the same intermediate type as observed in polygonal quan- 
tum billiards. Our central result is the evaluation of the spectral two-point 
correlation form factor at small argument, which in turn yields the asymp- 
totic level compressibility for macroscopic correlation lengths. 



1 Introduction 

The classification of quantum systems according to universal statistical proper- 
ties is one of the central objectives in the study of quantum chaos. It is generally 
believed that the spectral statistics of systems with chaotic classical limit are 
governed by random matrix ensembles, while systems with integrable classical 
dynamics follow the statistical properties of independent random variables from 
a Poisson process 0|2j|2]. Interestingly, certain billiards in rational polygons 
fall in neither of the two universality classes: the energy level correlations are 
conjectured to be of intermediate type El HI El ■ In particular, this means 
that the consecutive level spacing distribution P(s) exhibits level repulsion simi- 
lar to random matrix eigenvalues, but has an exponential tail as for independent 
random variables. Furthermore, the spectral form factor K^ij) is intermediate 
between and 1 in the limit r — 0. One standard example for a statistics of 
this type is the semi-Poisson distribution, for which P(s) — 4sexp(— 2s) and 
K 2 (t - 0) = 1/2. 

In this paper we present a one-parameter family of quantum maps whose 
spectral statistics are of a similar intermediate type as observed in polygonal 
billiards, cf. Figs . ^ and |21 The main result of our investigation is the evaluation 
of the spectral form factor K 2 {t) at small argument. It is based on a number- 
theoretic analysis which turns out to be considerably easier than the geometric 
approach required for billiards [7J- 
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Figure 1 : The consecutive level spacing distribution for Hilbert space dimension 
N = 7001 and a = 1/2 (left) and a = 2/3 (right). The Poisson distribution 
corresponds to P(s) = cxp(— s), semi- Poisson to P(s) — 4sexp(— 2s), and COE 
to the level spacing distributions of the circular orthogonal random matrix en- 
sembles. 
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Figure 2: The consecutive level spacing distribution for Hilbert space dimension 
N = 7001 and a = 3/5 (left) and a = 5/8 (right). 
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Figure 3: The consecutive level spacing distribution for a = (\/5 — l)/2 and 
Hilbert space dimensions N = 6997 (left) and N = 5867 (right) corresponding 
to values of e^/N ~ 32.1 and e^/N = 0.41, respectively. 
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Consider the following map of the two-torus T 2 = R 2 /Z 2 , 

•>■■■>•-■>•. (:) » <" 

where / is some 1-periodic function. The map is a concatenation $y = $o o pj 
of free motion $o and kick p f , 

The quantization of a torus map associates with it a unitary operator acting on 
the iV-dimensional Hilbert space of functions tp : 1m — * C with inner product 
= -/V -1 Eq=o i>*(Q)<l>(Q)- Here z ^ = z /^ z denotes the inte gers modulo 
N, and iV has the physical interpretation of an inverse Planck's constant. The 
quantum evolution operators f($o) and U(pf) corresponding to $o and pf, 
respectively, are defined by the matrix elements (cf. |H1 1101 ITT1 112) ) 

1 £z? 

(Q'\u@o)\Q) = nT,^ i- p2 + W Q)) > ( 3 ) 

p=o 

(Q / |^(p / )|Q) = <Q'|Q)e(-w(§)), (4) 

where V(q) is a periodic function defined by f(q) = —V'(q), and ejv(ac) — 
exp(2TTix/N) and e(x) — exp(27rix). Furthermore {/($/) = U($>o)U(pf) and 
thus 

(Q'|f/($/)|Q) = (Q'\U^ )\Q) e (-NV (§)) . (5) 

We are interested in the special case of the piecewise linear sawtooth potential 
V(q) = —a{q} for some real constant a, where {q} denotes the fractional part 
of q. In this case f(q) = a. The corresponding classical map $ a := $y is 
uniquely ergodic for irrational a (in particular, there are no periodic orbits) 
but not mixing. For rational a, the motion can be identified with an interval- 
exchange transformation. Note that in the momentum representation \P), with 
(Q\P) = N~ 1 / 2 eN(PQ), the operator C/($ Q ) has the representation 

(P'\U($ a )\P) = — e N (-P' 2 ) 1 ~ < Na ) (6) 

N ' l-e N (P- P' + Na) (J 

if Na £ Z and 

(P'|£/($ Q )|P) = e N (-P' 2 )(P'\P + Na) (7) 

otherwise. 

Since E/($ Q ) is unitary its eigenvalues are of the form exp(i#j) with eigen- 

phases < 6\ < 02 < • • • < 8m < 2ir; it is convenient to set 6q := 9m — 2ir. 



3 



The spacing distribution for consecutive levels is described by the probability 
density 



where the factor of N/2n ensures that s measures spacings on the scale of the 
mean level spacing 2tt/N. Figs. HTTl display the spacing distribution P(s) of 
the eigenphases of the matrix U(& a ) for both rational fFigs. lTl2^l and irrational 
(Fig. I3J values of a, and N a prime number. In the case of rational a = a/b 
one should avoid Hilbert space dimensions N divisible by b, since the matrix Q 
has highly singular statistics |181 119| . For irrational a, we find in Fig. [3] (left) 
a spacing distribution which resembles those of random matrices from the the 
Circular Orthogonal Ensemble (COE). COE statistics are normally expected 
for systems with chaotic classical limit and time-reversal symmetry El • In 
our case the time-reversal transformation which anti-commutes with U($ a ) is 
T" = U^q) 1 / 2 Tt/(<I>o) _1 / 2 , where T denotes the complex conjugation operator 
Tip := -0*. Localized spacing distributions of the type seen in Fig. |31 (right) occur 
when e, defined as the oriented distance of Not to the nearest integer, is of the 
order \/\fN. Such correlations arise as the perturbation of a rigid spectrum, 
and will be described in section [21 

The map <& a has in fact a further classical symmetry: it commutes with pi/ 2 
for any a. Thus, for N even, the corresponding operators U(§ a ) and U{p 1 / 2 ) 
commute and the eigenstates of U(Q a ) fall into two parity classes according to 
U{pi/2)tpj = <fij or U(pi/2)fj = —tfj, respectively. Our numerical experiments 
suggest that the level statistics for each subspectrum for even N are of the same 
type as those for odd N displayed in Figs. 11131 

A statistics which is more accessible from an analytical point of view is the 
two-point correlation density (which describes the distribution of all spacings) 







The Poisson summation formula applied to the m-sum yields 




(10) 



The spectral form factor is defined as the Fourier transform of i?2(s), 



K 2 (T) = -\Tr[U($ a r\ 



t = n/N. 



(11) 



In the following section we will calculate 




(12) 



= 1 
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where the limit is taken along suitable subsequences Ni, iVjj, oo of integers. 

To illustrate the relevance of this quantity, let us consider the counting function 
Af(L, £) for the number of eigenphases in the interval [£ — £ + t£] (mod 27r). 
The number variance is defined by 

^ {L) = h f*W(L,0-L} 2 <% (13) 

- # E E C -9 k + 2nrn)) L 2 (14) 

with = Lx(:e — y)x(y)dy — max{l — |x|,0}, where x is the indicator 

function of the interval [—5,5]. In view of |J3J, (|l(Jfl and 114|) . number variance 
and form factor are related by 

s2 W = Is E I Tr P(*«) B ] I 2 x (4 »Y ( 15 ) 

where x(y) = sin(7ry) / (ny) is the Fourier transform of \- It then follows from 
(|12[) and a standard probabilistic argument that, for macroscopic intervals of 
size L = £N (with I > fixed as N v — > 00), we have 

lim lim — i-i = if 2 (0) (16) 

since £^2 kjt0 x{k£) 2 -> Jx(y) 2 dy = 1, as ^ ->■ 0. The ratio T?{L)jL is called 
the ZeveZ compressibility 



2 Spectral form factor at small argument 

To calculate the value of K^{t) at small but non-zero values of r, we note that 
in the semiclassical limit N — ► 00, the trace TrJ7($ a ,) n is asymptotically equal 
to TtJ/($™), with n arbitrary but fixed. This may be seen by representing the 
respective traces as Gutzwiller-type periodic orbit sums. 

The advantage of the choice of <&„ over other piecewise linear maps is that 
there is an explicit formula for the nth iterate. (We have observed intermediate 
statistics also for the closely related "triangle maps" introduced by Casati and 
Prosen ^3]; the classical analysis of these maps is however considerably more 
involved.) Here, a short calculation shows that 

K =P(n-l)a/2°*0 °P(n+l)a/2- (17) 

The corresponding quantum evolution is therefore given by 

(Q'\U(K)\Q) = e (N (s=i) a {^}) (Q'\U(^)\Q)e (n (*±i) a {§}) , 

(18) 
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and so 



N-l 



N-l 



TrU(t>Z) = l£ e N {~nP 2 ) x ]T e(naQ) 



(19) 



p=o 



Q=Q 



where the first sum is a classical Gauss sum and the second a geometric sum. 
Let m := 2n/ gcd(2n, N) and M := N/gcd(2n,N), then 



N-l 



E ^(~nP 2 ) 



TV 
M 



M-l 



Eexp( m P 2 f; 



(20) 



which evaluates to N/VM if Mm is divisible by 2, and vanishes otherwise. The 
geometric sum is O(l) for irrational a and n ^ and hence K2(n/N) ~ for 
all bounded n in this case. For rational a = a/6, the geometric sum equals iV if 
n is divisible by 6 and is O(l) otherwise. Thus K2(n/N) ~ gcd(2n, TV) provided 
n is divisible by 6 and 2Nn/ gcd(2n, TV) 2 is divisible by 2; K2{n/N) ~ in all 
other cases. 

If we restrict ourselves to a subsequence of the values of N which are prime 
numbers then gcd(2n, N) = 1 for n < N, and the time averaged form factor is 



i " i 

K 2 (0) := lim lim - V K 2 (n'/N) = -. 

n^oo N^oo n ' 



N— >oo n 

iV prime n' — 1 



(21) 



These values of if 2(0) are consistent with those expected for intermediate statis- 
tics 0. The case b = 2 and N prime, for which ^(0) = 1/2, does however 
not agree with the Poisson statistics seen numerically in Fig. ^ (left), where 
^2(0) = 1. The solution to this apparent paradox is that P(s) displays correla- 
tions on the scale of the mean level spacing, whereas K 2 (0) involves correlations 
on much larger scales of order N. Similar discrepancies between a random 
matrix- like P(s) and a non-universal K 2 (t ~ 0) have been observed for non- 
arithmetic Hecke triangles JUE], compact hyperbolic triangles and tetrahedra 
[T?)] and cat maps coupled to a spin 1/2 [T7j . 

If N is twice a prime, i.e., N — 2R with R an odd prime, then gcd(2n, N) = 2 
for n < R, and 2Nn/A — Rn is divisible by 2 if and only if n is even. Hence 
only terms with n divisible by lcm(2, b) contribute, and so 



1 » 

K 2 (0) := lim lim - V K 2 (n'/N) = 



AT— >oo n 

N/2 prime 



lcm(2,6) 



(22) 



The situation is different for TV = pR, where p is a fixed odd prime, and R 
runs again over all odd primes. Now gcd(2n, N) = gcd(n, p) for n < R, and 
2Nn/ gcd(n, p) 2 is always divisible by 2. Since p is prime, gcd(n, p) = 1 if n is 
not divisible by p and gcd(n, p) — p if it is. In this case 



K 2 (Q) := lim lim - 

n^oo N^oo n 
N I p prime 



^W-i + H^- <23) 
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Figure 4: Left: Level compressibility E 2 (L)/L with L = 100 for rational a — 
a/b (x), versus K 2 (0) (©). From top to bottom: N = 7001 (Equation (|5T|l L 
TV = 6998 (Equation and TV = 6999 (Equation J23J). Right: nu merica l 

value of the form factor and theoretical plot l|24|) with ^(0) given by ^(0), 
for the same N and a = 1/4. 
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Fig. |21 (left) illustrates the asymptotic relation H16|) between the level com- 
pressibility and -^2(0) for large values of N: a prime (N = 7001), twice a prime 
(N = 6998) and three times a prime (N = 6999). Fig. [3 (right) compares the 
numerical form factor, obtained by diagonalizing the matrices £/($ a ), with the 
model |7| 

. . A 2 - 2A + 4tt 2 t 2 , . 

K ^ = A 2 + 4,V 2 ™ 



with A equal to 2/ (1— if 2 (0)). If N is divisible by b, the operator J7($ Q ) coincides 
with an alternative quantization Uo($ a ) := ^(^a/at) of <& Q proposed in |18|. 
where a is replaced by a rational approximation A/N so that |a — A/N\ < 1/N. 
The spectral statistics of Uo(§ a ) are well known to be highly singular and are 
not of intermediate type [IS]. It has been noted in |2U1[^ that Uo(Q a ) may be 
coupled to a spin 1/2 precession in such a way that intermediate statistics are 
seen numerically; the construction is analogous to the one for cat maps with 
spin 1/2 HZj. 




3 Localized level spacing distributions 

The above analysis yields trivially for irrational a 
1 » 

K 2 (0) := lim lim - V K 2 (ri IN) = (25) 

n' = l 

consistent with the COE statistics seen in Fig.|31(left). In contrast, Fig.|21(right) 
illustrates a class of statistics different from random matrix theory, which occur 
for subsequences of N, for which the quantity 

if{iVa}<!/2 
1 otherwise 

(the oriented distance of Na to the nearest integer) is at most of the order of 
l/y/~N. Note that if we take A/N to be the successive approximants in the 
continued fraction expansion of a irrational, we have e = 0(1/N). On the other 
hand, for rational a = a/b with N not divisible by 6, we have |e| > 1/6; the 
following considerations clearly to not apply in the latter case, where we may 
expect to see generic intermediate statistics. 

We shall now explain the localized level correlations observed for e = 0(1/ yN) 
by means of classical perturbation theory. The eigenphases 9^ , . . . , 9$ € 
[0, 27r) of Uo($ a ) and the corresponding orthonormal basis of eigenstates tp^ , . . . , 
are known explicitly, cf. |18) . Prop. 5.1. Since 

(Q'\U($ a )\Q) = (Q'\U (<f a )\Q) e N (eQ), (27) 

the Born expansion of the eigenphases 9j of C/($ Q ) is 9j = 9^ + e9^p + 0(e 2 ) 
with the first order correction given by 

flW _„,o„/,„(0)iai,„(0) 
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2M U '\A\^>), (28) 
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s s 

Figure 5: Integrated level spacing distribution for U(<f) (solid line) and inte- 
grated distribution of G(<p) for <f> uniformly distributed random variable (dashed 
line) for N = 431 (left) and N = 5867 (right). 



where A(</>) = {</>} — 1/2 is the sawtooth function. The term tt is irrelevant for 
the spacing distribution since it is independent of j. As to the second term, 
quantum unique ergodicity of Uo($ a ), proved in |18j . implies that in the limit 
N — > oo we have (<Pj | A\ifj) — > A(</>) d<f> = for all j. The eigenstates 

have a particularly simple form for N a prime number |18j . which we will 
assume in the following. In this case it can be shown that the level spacing 
distribution P(s) for the 9j is asymptotically given by the distribution of 

G{4>) = 1 + eVNj2 l ~ eN 2~£~ lk) e{H + M (29) 

fcGZ 

where (f> is a uniformly distributed random variable in [0,1), A~ x the inverse 
of A modulo N and (3k, n some explicitly known phase factor (see Appendix). 
The variance of the above distribution is e 2 A~ 1 (l — A^ 1 /N), with the choice 
A^ 1 e [0, N — 1]. Fig. [S] compares the numerical computation of the integrated 
level-spacing distribution J(s) = P(s')ds' for the eigenvalues 9j of U($ a ) 
(solid line) to the distribution of the random variable G((f>) (dashed line). 
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Appendix 



Let us define f3k,N by 

/ A ~ lk Ab N \ 1 N ~ 1 

(30) 

It can be shown that the variable £j := V~N(cpj \ A\ip^) is equal to £j — 
h(A~ 1 j/N), where h is the 1-periodic function 

h{x) = Y;&-ke{kx + p k , N ), (31) 
feez 

where Ao = 0, Afc = i/2nk (k ^ 0) are the Fourier coefficients of A(0). For iV 
prime, the spectrum of Uo(& a ) is totally rigid [T^], that is 6^ — 2nj /N+Cn for 
j = 1, . . . , N, where Cat is some overall shift. The corresponding level spacing 
distribution is hence P'°'(s) = 5(s — 1). For e small enough, the perturbation 
does not change the ordering of the levels. The level spacing distribution P(s) 
for the 9j is then given by the distribution of the (N/2Tr)(9j — 0j-i) ~ 1 + 
ev^V(^j — which is equal to G(A~ 1 j/N), where G is the function defined 

b ^ _ 

G(cj>) = 1 + eVivE A-fc e(/3 fe ,iv) [l - e N (-A- l k)} e(kcf>). (32) 

For random j 6 {l,...,iV} and Af large the distribution of P(s) is asymp- 
totically given by the distribution of G(0) where <j) is now a uniformly dis- 
tributed random variable in [0,1). The variance of the above distribution is 
(2eVN) 2 J2 ke z l A >=l 2 sin2 (^A~ l k/N) = e 2 A- l {l-A- 1 /N), provided we choose 
A- 1 £ [0,7V- 1]. 
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